Frequency spectrum of the Casimir force: interpretation and a 
paradox 
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Abstract. - The frequency spectrum of the Casimir force between two plates separated by vacuum 
as it appears in the Lifshitz formalism is reexamined and generalised as compared to previous 
works to allow for imperfectly reflecting plates. As previously reported by Ford [Phys. Rev. A 
48 (1993) 2962], the highly oscillatory nature of the frequency dependence of the Casimir force 
points to possibilities for very large and indeed negative Casimir forces if the frequency-dependent 
dielectric response, £(tj), of the materials could be tuned. A paradox occurs, however, because an 
alternative calculation of the effect of a perturbation of £(a;) involving a Wick rotation to imaginary 
frequencies indicate only very modest effects. A recent experiment appears to convincingly rule 
out the reality of Ford's optimistic predictions, although given the enormous technological promise 
of such frequency effects, further theoretical and experimental study is called for. 



In an interesting paper [1] , Ford analysed the frequency 
spectrum of the Casimir pressure as it appears when read 
directly out of Lifshitz' celebrated formula [2]. His cal- 
culations extended a previous study of the spectrum of 
the Casimir effect for a massless scalar field [3] and subse- 
quent analysis of the electromagnetic vacuum stress tensor 
by Hacyan et al. [4]. In [1], the classical Casimir set-up 
is considered, where two perfectly reflecting metal plates 
of infinite transverse size are separated by a vacuum-filled 
gap of width a. For this system, the pressure between the 
plates was found by Casimir [9] to be 
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Ford's puzzling finding was that if the pressure is expressed 
as an integral over all frequencies of the zero-point oscil- 
lations of the electromagnetic field in the cavity, the inte- 
grand is wildly oscillating and discontinuous as a function 
of frequency and the integral a sum of almost exactly can- 
celling positive and negative contributions, each of which 
far larger in magnitude than the measurable pressure it- 
self. By a suitable cutoff procedure, however, the integral 
is calculable and the result correct. Similar considerations 
were later performed for a sphere and plate set-up [5, 6] 
and the electromagnetic stress tensor in a cavity [7]. An 



extension of Ford's work on two ideally conducting plates 
was recently presented by Lang [8]. 

The unruly behaviour of the Casimir force as a func- 
tion of real frequencies was recently treated for numeri- 
cal purposes [10] and the same oscillatory behaviour was 
found. As a consequence these latter authors like most be- 
fore performed the Wick rotation to imaginary time (and 
imaginary frequencies), which is legal when the permit- 
tivity is assumed causal. As expressed for imaginary fre- 
quencies, the Lifshitz expression is much more well be- 
haved and rather than complex and strongly oscillating, 
the frequency integrand is real, nicely peaked and expo- 
nentially decreasing at high imaginary frequency. The im- 
portance of good optical data for the precise calculation of 
Casimir forces has been emphasised in a number of recent 
efforts [11-14], but these have all employed a Wick rotated 
formalism. 

Ford suggested that if the frequency response of the 
plate materials could be tuned, for example if a mate- 
rial could be found which is transparent for all but a 
small band of frequencies in which it was a good reflec- 
tor, Casimir forces much larger than that between perfect 
conductors could be observed and by changing the reflec- 
tion band the force could be changed from attractive to 
repulsive. Despite the potentially enormous technological 
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potential of such tuning for applications in microengineer- 
ing and nanotechnology, the issue of the physical interpre- 
tation of the frequency spectrum of the Casimir force has 
remained largely unaddressed. 

A relevant experiment was recently performed by lan- 
nuzzi, Lisanti and Capasso at Harvard [15] in which the 
Casimir force was measured in the same configuration with 
a good and poor reflector respectively. The material used 
was a so-called hydrogen-switchable mirror which can be 
switched from mirror to transparent at optical frequencies 
by introducing hydrogen. At frequencies in the IR and 
UV parts of the spectrum, reflection is presumed by the 
authors of [15] to be approximately unchanged. This is 
the inverse of the situation suggested by Ford, and based 
on [1] large effects should be expected. lannuzzi et al. 
observed no change of the force in the two cases, however. 

The paradox is theoretical as well as experimental: As- 
sume that the permittivity of one or both of the plates in 
a standard two-plate set-up [9] is changed in a band of fre- 
quencies. The effect on the force of this perturbation may, 
one may think, be calculated in two different ways. Either 
Lifshitz' expression for the force is integrated over the rel- 
evant band of real frequencies and the difference taken. 
Alternatively, the new, perturbed permittivity is rotated 
to imaginary frequencies using the Kramers-Kronig rela- 
tions and the force calculated as an integral (assuming 
zero temperature) over imaginary Matsubara frequencies. 
The numerical results obtained in the two ways, however, 
appear to differ greatly. 

In the following. Ford's calculations are generalised to 
imperfectly reflecting plates represented by an effective re- 
flection coefficient. While the same oscillatory behaviour 
is found, the integrand is no longer discontinuous. In the 
limit of perfect reflection, Casimir's result is once more 
obtained. Thereafter the paradox is elaborated and com- 
pared to the recent experiment by lannuzzi et al. 

Casimir force and frequency spectrum with con- 
stant reflection coefflcients. — The expression for the 
Casimir pressure was given by Lifshitz as an integral over 
(real) frequencies and the variable p by [2] 



show that the former part of the integral covers all modes 
which propagate in vacuum, that is when lu > |kj^|c, and 
the latter covers all evanescent modes lo < |kx|c which 
vanish exponentially away from the surfaces. 
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Fig. 1: The integration contour in (2). 

In the first part of the analysis we shall consider complex 
reflection coefficients which are constant with respect to 
to and p, but in principle dependent on the separation a. 
We will assume zero temperature throughout so that the 
coth function in (2) is set to 1. 

To denote the frequency dependence we use the notation 
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Let us assume now that the reflection coefficient is a 
constant quantity with respect to lo and p which we allow 
to be complex for generality. 

Consider first the propagating part of (2). The p inte- 
gral may now be evaluated explicitly. Consider one of the 
modes and suppress the index a for now to find 
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where 7 = and ^ = 2uja/c. We now introduce the 
polylogarithmic function whose mth order is defined 
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Then we may write the pressure contribution from prop- 
agating waves (PW) P^^ as 



exp(2ipa;a/c) 



(2) 



where the variable p relates to the transverse wave-vector 
kj_ as 

p = — A/k^ — uj'^/c^ 

and the contour of integration C, shown in figure 1, runs 
from 1 to along the real axis and thence along the imag- 
inary axis to +ioo. We have assumed both plates equal in 
(2) but all results in the following may be generalised to 
different media 1 and 2 by letting — > rir2- It is easy to 
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When r„ are real, there is no evanescent contribution to 
the pressure as we will see. We have plotted the integrand 
(3) as a function of ^ for various real values of the effective 
reflection coefficient r in figure 2 (assumed equal for TE 
and TM for simplicity). In the limit r = 1 one obtains 
the discontinuous behaviour reported by Ford. Equation 
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Fig. 3: The Casimir pressure for non-unitary real reflection 
Fig. 2: Plot of (3) normalised by as a function of a;a/c for coefiicient (assumed equal for TE and TM) plotted relative to 
different real values of r. This figure generalises figure 2 of [1]. Casimir's result (1). 



(3) thus generalises Ford's calculation; indeed figure 2 is 
formatted so as to ease comparison with figure 2 of [1] . 

Wc can go on to calculate the force for constant reflec- 
tion coefficients. The analysis so far only required to be 
independent of p; one could define some p averaged r„{u)) 
and plot it the way we have in figure 2. By assuming r 
independent of lo as well, however, we are able to evalu- 
ate the w integral explicitly. We substitue w by ^ in the 
integral and write Ua- = exp(i^) for short. 

For a ^ independent one may show from the recursion 
property [16] 



remains: 
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that the following relations are fulfilled: 



Jo 

Jo 



Re 



The integrand as given by (3) is a wildly oscillating func- 
tion of ^ as figure 2 indicates. For dcfinitcncss, let r 
rexp(— (5^) where 5 is a small real quantity which we will 
take to zero in the end. The same procedure was followed 
by Ford [1] and reflects the physical fact that all materials 
become transparent at very high frequencies. By moder- 
ately lengthy but straightforward partial integrations of 
(3) with respect to ^ one may then obtain for constant 
and real reflection coefficients 
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With the normalising factor exp(— 5^) in all reflection co- 
efficients the oo limit is zero as physically expected. In the 



limit ^ = 0, Ua- 



and only the term with no factor ^ 
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Secondly we consider the part of the pressure from 
evanescent waves (EW). The same way as above we eval- 
uate the p integral thus: 
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where p = iq. Clearly if is real as in a non-dissipative 
medium, there is no evanescent contribution. Moreover, 
one notices that that the evanescent contribution is ex- 
actly cancelled by the last term of (3) above. Using the 
same normalisation as before the evanescent part accumu- 
lates to zero. Note that we could have let the normalisa- 
tion exp(— 5^) pertain to the exponential factor exp(i^) in 
the integration (4) since the remaining divergent term is 
exactly cancelled by evanescent contributions. 

Hence the final result is obtained when is constant: 
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In the limit of perfect reflection, ^ 1, the summation 
over a gives a factor 2, and with Li4(l) = ^(4) = 7r^/90 
we get Casimir's result (1). P{a,r) is plotted in figure 3 
relative to the case for ideal reflection. 

In exactly the same manner we could find the Casimir 
free energy for real and constant reflection coefficients to 
be given as 
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Again Casimir's result is obtained in the limit = 1. 
Some caution should be excerted here since function Jji4{z) 
has a branch cut along the real axis from 2; = 1 to oo across 
which its imaginary part is discontinuous. Its real part, 
however, is continuous everywhere [16]. 
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Fig. 4: rcr{a) for gold as calculated from calculated force using 
(8). 



Physical significance: naive approach. Phys- 
ically, reflection coefficients are far from constant with 
respect to frequency and transverse momentum, but the 
above analysis using constant reflection coefficients can be 
thought of as an averaging process in which the constant 
To- is defined as the value which, when replacing the real- 
istic ra-{oJ,p) does not change the value of the integral (2). 
Such an interpretation requires to depend on a as well. 

Using the calculated Casimir pressure between gold 
plates where realistic data from [17] are employed (the 
calculation was previously presented in [19]), we can cal- 
culate the effective reflection coefficients To- (a) by 
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where Pnum are the calculated data, and Li~^ is the in- 
verse polylogarithm of order m. By use of the series re- 
version functionality of analytic software the inverse poly- 
logarithm is simple to calculate numerically. We plot the 
effective reflection coefficient in figure 4. 

Ford [1] (and later Lang [8]) appears to argue that the 
frequency spectrum shown in figure 2 resembles the real 
spectrum, with the modification that the oscillations will 
be dampened and finally vanish at high frequencies. He 
concludes from this that since the oscillations themselves 
are much larger than the final force (which is the remain- 
der of large fluctuations which cancel each other almost 
exactly), much larger and even repulsive Casimir forces 
could be obtained if media could be found which were 
good reflectors only in a range of frequencies. 

Such a simple physical interpretation of the frequency 
spectrum is not unproblematic however. For comparison 
with a more physically realistic system we have plotted 
the frequency spectrum obtained by plugging tabulated 
optical data for gold from Palik's book [17] directly into 
(2) and integrating over transverse momentum numeri- 
cally (the separation is lOOnm; the integrand is now a 
function of lo and a individually, not only their product), 
as well as the same integrand as obtained when using the 
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Fig. 5: The frequency spectrum of the Casimir pressure cal- 
culated using real data for gold and the plasma model respec- 
tively. We have used a — lOOnm. The below graph is a zoom 
of the above. 



simple plasma model e(a;) = 1 — o^p/tj^ using plasma fre- 
quency for gold uip = 9eV. Although the force predicted 
by either model at zero temperature are found to be very 
similar in magnitude when calculated by way of Wick ro- 
tation^ , the frequency spectrum clearly differs greatly. 

While a useful generalisation, the formalism of mean 
reflection coefficients should be used with care since its 
substitution renders the integrand of (4) void of physical 
meaning other than giving the correct value per definition 
after integration: u simply becomes a dummy variable. 
If Ford's result can be seen as a limiting case of the same 
procedure, this would indicate that the spectrum depicted 
in figure 2, while interesting may not represent physics. 

Nonetheless, the key feature of the integrands of figures 
2 and 5, the presence of large oscillations, is a hallmark 
of all these graphs. These fluctuations should be physi- 
cally observable as noted by several authors. As we will 
see, however, indications are that while experimental con- 
firmation of the Casimir force as calculated using Wick 
rotation is plentiful (see [18] for a review), calculation by 
straightforward integration over real frequencies seems at 



^Disregarding the disputed behaviour of the TE mode near zero 
frequency [19]. 
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odds with a recent experiment. This is paradixical since 
the two methods are typically presented as equivalent. 

The paradox and an experiment. Assume that 
the permittivity of two metallic planes in a Casimir cavity 
is perturbed in such a way that it is made transparent in 
a band of frequencies but is still a good reflector outside 
this band (with reflectivity dying ofl^ at high frequencies). 
The effect on the force from this perturbation may be cal- 
culated in two different ways and the results are different. 

Exactly such a situation was probed experimentally by 
lannuzzi et al. [15]. The force between two metallic plates 
(in reality a sphere and a plane) was measured, one of 
which a hydrogen-switchable mirror (HSM). A good mir- 
ror in its as-deposited state, the HSM becomes transpar- 
ent in the visible region once in a hydrogen rich atmo- 
sphere. According to the authors of [15] the switching of 
the mirror corresponded roughly to setting the material 
reflectivity to zero over a wavelength range 0.2 — 2.5 fim, 
corresponding to co between about coi = 7.5 • lO^^l/s and 
UJ2 = 9.4 • 10^^1/s. At a precision of 10-15% (measured 
roughly from figure 4 in [15]) the group were unable to 
detect any difference in the force with the mirror switched 
on and off respectively. 

We will estimate the effect of the transparency window 
using two different methods (an idealised version of such 
a material was considered in [8]). Assume first that the 
boundaries of the transparent window, between frequen- 
cies Wi and UJ2, are sharp so that for a complex permittivity 
e = e' + ie" = 1 + X, 

X(w) X('^)[l - 0{u; - ui)0{oj2 - w)] 

where 8 is the unit step function. Using, as in [15], the 
Drude model with data for gold, e{ui) = 1 — uip/{ui'^ + 
icou) with iVp = 9eV and u = 35meV, the change in the 

permittivity at imaginary frequencies is fond as e(iC) 
e{iC) — ^e(*C) by use of the Kramers Kronig relation 
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the use of which ensures that the perturbation obeys 
causality. Ae(i^) makes for a correction to e{iQ) on the 
level of l%-4% and a corresponding correction to the force 
which would be unobservable at the precision of the exper- 
iment [15]. The latter authors use precisely this argument 
to explain their negative result. 

Now let us calculate the correction by instead inserting 
a modified e(a;) into (2) where for a single interface 
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and instead of the step functions we model the transparent 
window using a function which allows smooth boundaries. 
Let x(w) -> x('^)'/>('^; A, s) with 
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Here the parameter A G [0, 1] is the relative reduction in 
the window and s determines the sharpness of the edges 
with s = 00 giving unit step behaviour. The function ^ is 
plotted in figure 6 for A = 1 and different values of s. 

We now calculate the change in the force by simply in- 
serting this modified e{w) into (2) and integrate over p 
and a sufficiently large frequency range, then taking the 
difference. We use the Drude and plasma model respec- 
tively with parameters for gold to model e (the data of [17] 
contain too few points to use without extensive extrapo- 
lation). In the experiment [15] only one of the plates was 
gold, yet the calculations are so rough that the difference 
does not matter here. The fact that in this estimation 
both refiection coefficients vanish in the transparency win- 
dow whilst in the experiment only one did so, makes for 
a slight overestimation of the effect for smooth window 
edges whilst it makes no difference when edges are sharp 
since the force only depends on the product of the reflec- 
tion coefficients of the two materials. Note furthermore 
that the erratic behaviour of the integrand makes the nu- 
merical accuracy of integration somewhat rough, and also 
that in the modelling of the transparency band no effort 
has been made to ensure causality is satisfied. Thus the 
estimate is accurate only to order of magnitude. 

The results for different values of A and s at separa- 
tion a = lOOnm are shown in figure 7. Note that the 
absolute pressure between parallel gold plates at T = is 
approximately 6 Pa (e.g. [19]) and that the difference be- 
tween using Drude and plasma models is negligible here. 
Clearly the corrections at A close to 1 are much too large, 
more than 10 times the force itself. The smoothness of 
the edges has no obvious effect. Only the plasma model 
with sharp boundaries and A < 0.5 could fall within the 
experimental accuracy of [15], yet it seems highly likely 
that the realitvely small corrections at low A and high 
s for the plasma model is due to chance, especially since 
the more realitic Drude model gives corrections which are 
enormous and also, counterintuitively, positive. Such large 
corrections would be in keeping with Ford's predictions, 
but seem clearly ruled out by experiment. 

Conclusions. — We have revisited the question of 

the frequency spectrum of the Casimir force, generalising 
Ford's result [1] to the case of subunitary reflection coeffi- 
cients. While this smooths out the frequency spectrum as 
given by the real frequency Lifshitz formula integrand, the 
integrand is still wildly oscillating. This result is used to 
calculate the Casimir force and free energy for constant, 
"effective" , reflection coefficients, a new result to the au- 
thor's knowledge. 



p-5 



S. A. EUingsen 



(p (q) 



1 




















8- 


,\ 
\\ 








/ - 
/.' 









6- 
4- 
2- 
0- 








■'i 




— s-500 

--■ 3 = 100 

— ■ 20 
s- 10 








// 
.-' / ( 

— -J 







12 3 4 

Qa/ c 



Fig. 6: •-piyi) for A = 1 with different s values. 

While the effective reflection coefficient may be a use- 
ful model, there is little reason to believe the resulting 
frequency spectrum to represent physics. A more realistic 
spectrum using physical models for the permittivity of ma- 
terials is still highly irregular and, taken at face value, in- 
dicates that large and even repulsive Casimir forces could 
be attainable by tuning the dielectric response of materials 
used. An alternative means of calculation, paradoxically, 
gives a different, pessimistic result, and the large effects 
seemingly implied by the wildly behaved frequency spec- 
trum will seem to be excluded by a recent experiment by 
lanuzzi et alia. 

With the possibility of technological applications of the 
Casimir force however, there is reason to strive for a bet- 
ter understanding of the physical interpretation of the fre- 
quency spectrum of the Casimir force as well as further 
experimental efforts to settle this issue. An experiment 
similar to [15] accompanied by a careful measurement of 
the dielectric response of the actual sample used over a 
large frequency region would be a straightforward possi- 
bility, and a more sensitive measurement of the force might 
also be able to measure the actual difference in pressure. 

It appears that the straightforward interpretation of the 
Casimir frequency spectrum as the integrand of the Lif- 
shitz force formula at real frequencies is not valid, yet 
given the modesty of the efforts presented herein further 
investigation is warranted. Furthermore, the paradox pre- 
sented herein pends a satisfactory resolution, hopefully to 
appear in the future. 

The author thanks Irina Pirozhenko for supplying the 
data from [17] in electronic format and acknowledges 
highly useful suggestions from an anonymous referee. 
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Fig. 7: Force difference with and without a transparent window 
using direct frequency calculation. Circles: Drude model, no 
circles: plasma model. 
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